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Abstract 

The probability that a randomly accelerated particle in two dimensions has 
not yet left a simply connected domain A after a time t decays as e~ Eot for 
long times. The same quantity Eq also determines the confinement free energy 
per unit length A/ = fcgT Eq of a semiflexible polymer in a narrow cylindrical 
pore with cross section A. From simulations of a randomly accelerated particle 
we estimate the universal amplitude of A/ for both circular and rectangular 
cross sections. 
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Consider a long semiflexible polymer with persistence length P fluctuating in a cylindrical 
pore with diameter D. In the narrow-pore limit D « P the free energy of confinement per 
unit length Af is given by 

^f = A Op^rr (1) 



This follows from simple scaling or dimensional arguments [T]-[3[, such as that given below. 
Similarly, for a pore with a rectangular cross section || with edges Li,L 2 << P 

A ^ = A n ~j=k/3 ( + ) • ( 2 ) 

The dimensionless constants Aq, A^ in equations ([I]) and (H) are universal numbers, 
independent of both macroscopic and microscopic properties of the polymer chain. From 
computer simulations Dijkstra et al. estimated 

A Q = 2.46 ±0.07. (3) 

Solving an integral equation numerically that arises in an exact analytical approach, 
Burkhardt obtained || 

A u = 1.1036. (4) 

In the Monte Carlo simulations of Dijkstra et al. f2|, polymer configurations consistent 
with the Boltzmann distribution were generated with a Metropolis algorithm incorporating 
detailed balance. In this Letter we estimate Aq, A^ with a precision of about 1% using a 
Langevin dynamics approach, which is simple and efficient. Instead of a confined semiflexible 
polymer, we simulate a Newtonian particle which is randomly accelerated by Gaussian white 
noise in two dimensions. In the narrow-pore or tight- confinement limit these two systems 
have equivalent statistical properties, as discussed in and reviewed below. The basic idea, 
illustrated in figure 1, is that each of the possible world lines of a particle, which is randomly 
accelerated in two dimensions and remains in a domain A for a time t, corresponds to an 
allowed configuration of a tightly-confined semiflexible polymer in a cylindrical pore or tube 
with cross section A. 



We now describe the correspondence in more detail. As in ||, polymer configurations 
are specified in terms of Cartesian coordinates (x, t) = (xi,X2,t). The t axis is parallel to 
the axis of the cylindrical pore, as shown in figure 1. In the narrow-pore limit D << P or 
Li,L2 << P, configurations with overhangs are negligible, i.e. x is a single-valued function 
of t. The partition function is given by the path integral |3| 

Z(x, u; x , u ; t) = D 2 x exp [ - -P^ dt I — J (5) 

where x and u = dit/dt denote the displacement and slope of the polymer at t, and xo and 
Uo the same quantities at t — 0. A hard- wall confining potential is assumed, and the x 
integration is limited to the domain A. 

The path integral implies the the partial differential equation 

J^ + u- V x - ^V u 2 ] Z(x,u;x u ;t) = (6) 

to be solved with the initial condition 

Z(x, u; x , u ; 0) = 5(x - x )5(u - u ). (7) 

Since the polymer is confined to the interior of the pore and configurations with a disconti- 
nuity in slope cost an infinite energy, Z(x, u; x , u ; t) vanishes for u • n > as x approaches 
a hard wall. Here n is a vector normal to the wall and directed toward the interior of the 
pore. 

As in [^] we consider exponentially decaying solutions of equation (||) with the form 
u) exp(— Et). The eigenf unctions ip and eigenvalues E n satisfy 

u ■ V x - — V u 2 - E n ] ^ n (x, u) = (8) 

where ^(x, u) also vanishes for u ■ n > as x approaches the pore wall. In the long-polymer 
limit the partition function and the confinement free energy per unit length are given by 

Z(x, u; x , u ; t) ~ const x ipofe., u)-^ ( x o, — uo)e~ Eot t — > oo (9) 
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j^ = W) (10) 

for a pore with a circular cross section of diameter D, where E is the smallest of the 
eigenvalues E n . 

The P and D dependence in equation (JJ) may be derived from fllPf) by making the scale 
change x = ax', u = /3u' in where a and (3 are arbitrary positive constants. This leads 
to the scaling property 

E n (P, D) = a- 1 /3S n (a- 1 /3 3 P, a^D) (11) 
of the eigenvaues. Setting a = D and /3 = P~ 1 / 3 D 1 ^ 3 in equation ( |TTD and substituting the 



result in (jT0|), we obtain 

A/ P (M) 



(12) 



k B T PV3D2/3 

This expression is entirely consistent with equation (H) and implies 

A Q = E (1,1). (13) 

Similarly, for a pore with a rectangular cross section 

2A n = £ (l,l,l) (14) 

where E (P, L ly L 2 ) denotes the smallest eigenvalue in equation (|p for a rectangular domain 
with edges L\, L 2 . Equations (|I3[) and (|T^) allow us to determine the universal amplitudes 



Aq and A^ from calculations with P — D — L\ — L 2 — 1. 

For the rectangular domain equation (||) has separable eigenfunctions ip m , n (x,u) 
m (x 1 , Ml )0 n (x 2 ,M 2 ), with eigenvalues P m , n = ^ ldim )(P, L x ) + E^ dim \P,L 2 ), where 

9 1 <9 2 



-P,( n ldim )(P,L) m ( a; ,ii)=O. (15) 



<9x 2P 9m 2 

Here < x < L, and 4>(x,u) vanishes for x = 0, u > and x = L,u < 0. Thus equation 
(|14]) may be rewritten as 



^ D = ^ ldim) (l,l) (16) 

in terms of the smallest eigenvalue of ( |I5| ) for P — L = 1. The earlier numerical result for 
A^, noted in equation (||), was obtained in by converting (|15|) to an integral equation, 
determining the smallest eigenvalue numerically, and substituting the result in (p~6|) . 

Now consider a particle which is randomly accelerated in the <i-dimensional space 
(xi, . . . Xd) by Gaussian white noise with zero mean according to 

^ = m (mm) = p~%Kt - f). (i7) 

The probability density in phase space (x, u) that the particle remains in the domain A for a 
time t while the position and velocity evolve from (xo, uo) to (x, u) satisfies a Fokker-Planck 
equation || with exactly the same form (||), initial condition (]?]), and boundary condition 
at the boundary of A. Thus the probability density of the randomly accelerated particle 
equals the partition function of a semiflexible polymer in a pore with cross section A and 
decays as exp(— E t) for long times, as in equation (|]). 

To estimate A^ using equation (|1^), we simulated a randomly accelerated particle in one 
spatial dimension. From simulations in two dimensions we obtained a second estimate of A^ 
based on ( p!4[ ) and a prediction for Aq from flT3|). The simulation routine, similar to that in 
||, will now be described briefly. 

In an unbounded d- dimensional space the exact solution of the Fokker-Planck equation 
(P) with initial condition (0) and with P = 1 is given by [p,10[ 

Z free (x, u; x , u ; t) = (^^j ex P (x - x - u t) • (x - x - ut) + i(u - u )H 2 j 

(18) 



Trajectories with the probability distribution Zf rco (x n+ i, u n+ i; x n , u n ; A n+i ) given by ([18] 
are generated using the algorithm 



x n+1 = x n + u n A n+ i + e n £?n+i (3 1/2 s n +i + r„+i) (19) 
u n+1 = u„ + e n 2Ai^ r n+1 (20) 



where x n and u n are the position and velocity of the particle at time t n , and A n+ i = t n+ \ —t n . 
The quantity e n is a unit vector that points either along the positive X\ axis , the positive 
22 axis, . . ., or the positive Xd axis with equal probability, and r n and s n are independent 
Gaussian random numbers satisfying 

(r n > = ( Sn >=0, (r£> = <<£> = I- ( 21 ) 

In the absence of boundaries there is no time-step error in the above algorithm,i.e., the 
A n may be chosen arbitrarily. Close to the boundaries small time steps are needed. As in 
we performed our simulations with 

A n+1 = 10~ 5 + 10- 1 D n (22) 

where D n is the distance from the particle to the closest point of the domain boundary at 
time t n . This time step fulfills a reliability criterion discussed in 0. 

Some sample simulational results are shown in figure 2. The quantity Q(t) is the prob- 
ability that a particle with a random initial position in a one or two-dimensional domain 
and with initial velocity zero, which is randomly accelerated according to equation ( |TTD with 
P = 1, has not yet left the domain after a time t. The curves labelled 'circle,' 'square,' and 
'interval' refer to a circular domain of diameter D = 1, a square with edges Li = L 2 = 1, and 
a one- dimensional interval of length L — 1. Each of curves is based on 10,000 independent 
trajectories. 

Since Q{t) = J d d x J d d u J d d x Z(x, u; x , 0; t), it decays as e~ Eot for long times, with 
the same decay constant E as Z(x, u; x , u ; t) in equation (||). We determined E for 
the circular, square, and one-dimensional domains by fitting Q(t) for long times with an 
exponential function. A surprising result, shown in figure 2, is that the curves for the 
circular and square domains practically coincide when plotted versus E t instead of t. 

We also estimated E for circular, square, and one-dimensional domains from trajecto- 
ries that all begin at the the center of the domain with initial velocity zero. The results 
are consistent with the results for random initial positions but have a somewhat greater 
statistical uncertainty. 
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Combining our best estimates of E with equations (|13|), (0), and fll6|), we obtain 

A Q = 2.375 ± 0.013 (23) 
A u = 1.108 ± 0.013 (24) 

The result (E3I) for Aq is somewhat lower than the estimate (^) of Dijkstra et al. |2j 
and has a smaller statistical uncertainty. The result (|24|) for A^ in equation is consistent 
with the value (f|) obtained by Burkhardt from the numerical solution of an exact integral 
equation ||. Bundschuh [5 L] has also confirmed (£|) to within a few percent with a numerical 



transfer matrix approach |L2| for semiflexible polymers 



In summary we have found a simple and efficient simulational procedure for calculating 
the confinement free energy of a semiflexible polymer in a narrow cylindrical pore with cross 
section A. We use the equivalent statistical properties of a Newtonian particle which is 
randomly accelerated by Gaussian white noise in two dimensions. The probability that the 
particle has not yet left a domain A in a time t decays as e~ E ° l . We determine Eq from 
our simulations and then interpret it, following equation fliTf), in terms of the polymer free 
energy. We emphasize that the equivalence between the statistics of the polymer and the 
randomly accelerated particle is asymptotically exact in the limit in which the pore diameter 
is much smaller than the polymer persistence length. 

TWB appreciates the hospitality of the Institut Laue-Langevin, where this work was 
begun, and thanks Ralf Bundschuh for helpful correspondence. 
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FIGURES 




FIG. 1. The curve may be interpreted as a tightly-confined semiflexible polymer in a cylindrical 
pore with cross section A or as the world line of a particle which is randomly accelerated in two 
dimensions and remains in a domain A for a time t. 
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FIG. 2. Probability Q(t) that a particle with random initial position in a one or two-dimensional 
domain and initial velocity zero, which is randomly accelerated according to equation (17) with 
P = 1, has not yet left the domain after a time t. The curves labelled 'circle' (solid line), 
'square' (dashed line), and 'interval' (solid line) correspond to a circular domain with diameter 
1, a square domain with edge 1, and a one-dimensional interval with length 1. For the three curves 
Eq = 2.375, 2.199, 1.108, respectively. The dot-dashed line represents a pure exponential decay 
e~ Eot . 
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